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1.  INTRODUCTION 


This  report  discusses  the  theory,  applications  and  performance  of  various 
stages  of  analysis  and  data  processing  involved  in  the  DAASM  (Doppler  arrival 
angle  spectral  measurements)  Project.  The  analysis  can  be  divided  into  three 
stages,  namely: 

1)  Power  spectral  estimation 

2)  Coherence  analysis 

3)  Doppler  frequency  - angle  of  arrival  estimation  (DAASM  maps) 

The  power  spectral,  coherence  and  DAASM  maps  are  presented  in  the  form  of  digi- 
coder  printouts.  The  printouts  display  integers  0 to  15  in  a gradation  of 
increasing  intensity. 

Much  of  the  analysis  pertaining  to  the  powjr  spectra  and  coherence  has  been 
dealt  with  in  detail  in  previous  AFCRL  reports.  Only  brief  and  necessary 
material  are  discussed  here.  The  main  portion  of  this  report  deals  with  the 
adaptive  data  processing  aspect  of  the  DAASM  maps.  Several  test  functions  have 
been  devised  to  simulate  sources  at  different  Doppler  frequencies  and  different 
arrival  angles  and  with  various  levels  of  isotropic  white  noise.  The  simulation 
ir  performed  on  linear  arrays  of  six  and  twelve  antenna  elements  corresponding 
to  most  often  used  configurations  for  measuring  and  gathering  real  data. 

The  work  described  here  has  evolved  out  of  a need  to  extract  the  Doppler 
frequencies  and  the  direction  of  arrival  of  high  frequency  radio  signals  propa- 
gating through  the  ionosphere  and  received  on  a linear  array  of  antennas.  The 
concept  of  a homogeneous  random  electromagnetic  field  in  the  sense  of  Yaglora[5] 
propagating  across  an  array  of  sensors  has  been  employed.  If  the  random  field 
can  be  considered  as  a wide-sense  stationary  process  it  is  possible  to  apply  a 
multi-dimensional  Fourier  transformation  and  obtain  the  so-called  Frequency  - 
Wavenumber  Power  Density  Spectrum  (hWPD).  The  frequency  dimension  resolves 
the  Doppler  or  temporal  components  while  the  wavenumber  dimension  yields  the 
direction  of  arrival  or  the  spatial  components  of  vhe  array  data. 

The  frequency-wavenumber  approach  has  been  developed  and  widely  used  in 
the  processing  of  seismic  signals  by  Capon  et  al.  [1].  All  spectral  analysis 
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techniques  utilize  some  or  other  form  of  data  windows  either  in  the  original 
untransformed  domain  or  in  the  transformed  domain.  Data  windows  are  necessary 
to  reduce  the  statistical  variability  of  the  spectral  estimate  and  also  to 
suppress  sidelobe  leakage.  As  a consequence  the  choice  of  data  windows  has  to 
compromise  between  conflicting  requirements  of  spectral  resolution  and  statisti- 
cal stability.  Besides,  the  data  windows  in  usual  spectral  analysis  are  of  a 
fixed  nature  in  that  they  do  not  depend  on  the  data  or  noise  characteristics. 

The  data  window  is  a taper  function  applied  to  the  estimated  correlation  func- 
tion and  the  product  is  transformed  to  get  the  spectrum.  Consequently,  the 
estimated  spectrum  is  smeared  by  the  convolution  of  the  window  function  and  the 
true  spectrum.  The  selection  of  a suitable  data  window  becomes  very  critical 
when  only  limited  data  are  available  or  cannot  be  processed  because  of  station- 
arity  considerations  and  adequate  resolution  has  to  be  maintained. 

Spectral  techniques  which  are  data-adaptive  and  offer  higher  resolution 
than  usual  methods  have  been  increasingly  used  in  seismology  as  mentioned  above. 
The  maximum  likelihood  method  of  estimation  of  the  frequency-waventunber  spectrian 
which  is  employed  in  this  study  belongs  to  the  data-adaptive  high  resolution 
class  of  spectral  estimators.  The  method  as  applied  to  multi-channel  streair.s 
of  complex  valued  data  points  from  a two-dimensional  array  of  antennas  in  the 
Doppler-arrival -angle-spectral  measurements  (DAASM)  experiment  has  been  discussed 
and  derived  in  detail  in  a;'  AFCRL  report  by  (Pfister,  Sales  and  V^vad)  [2] . 

There,  the  method  is  rendered  data-adaptive  in  both  the  temporal  as  well  as  the 
spatial  domains.  But  it  suffers  from  the  disadvantage  that  the  core  and  time 
requirements  for  computer  processing  become  unrealistic  as  the  number  of  antennas 
and  time  data  points  are  increased.  It  is  probably  feasible  on  such  giant 
parallel  computer  systems  as  the  Iliac  IV. 

Moreover,  the  improvement  observed  in  the  frequency  dimension  does  not 
justify  the  additional  complexity  contrasted  with  simpler  and  much  faster  methods 
utilizing  the  Fast  Fourier  Transform.  Tlierefore,  the  technique  as  employed  here 
is  adaptive  only  in  the  spatial  domain,  since  it  is  conceivable  that  the  number 
of  available  spatial  data  points  is  normally  far  less  than  the  number  of  tiir,e 
domain  data  points. 
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2.  DESCRIPTION  OF  THE  EXPERIMENT 


2.1  Overall  Project 

The  DAASM  system  consists  of  a twenty-element  array,  eleven  of  which  are 
to  form  a minimum  redundancy  (MR)  configuration  as  described  by  Richards  [3]  in 
a previous  AFC'RL  report.  This  neans  that  by  proper  selection  of  antennas  it  is 
possible  to  obtain  spacings  equal  to  successive  multiples  of  an  available  minimum 
spacing.  It  can  also  be  used  in  quasi-minimum  redundancy  configurations  where 
the  available  spacings  are  near  successive  multiples.  The  advantage  of  the  MR 
array  is  that  it  is  equivalent  to  a large  filled  array  with  the  same  aperture 
and  the  same  basic  spacing.  The  signals  incident  on  the  antennas  are  sampled 
in  time  and  recorded  in  the  form  of  complex-valued  digital  data  after  proper  IF 
conversion  of  the  carrier  frequency.  The  pulses  at  chosen  carrier  frequencies 
originate  either  /rom  a flying  aircraft  or  from  a ground-based  transmitter. 

The  object  of  subsequent  processing  of  the  recorded  data  is  to  extract  the 
spatial  and  temporal  characteristics  of  the  signals.  To  describe  tliese  charac- 
teristics it  is  necessary  to  introduce  the  concept  of  a frequency-wavenumber 
power  density  (FWPD)  spectrum.  In  the  case  of  acoustic  and  seismic  signals, 
wavenumber  is  a quantity  which  describes  waves  in  terms  of  their  velocities, 
that  is,  speed  and  direction,  are  both  variable.  But  for  electromagnetic  waves 
wavenumber  is  only  dependent  on  direction  since  the  speed  is  constant  in  a given 
homogeneous  medium. 

If  the  time  samples  are  understood  to  constitute  one  dimension  and  the 
spatial  samples  to  be  the  second  dimension  then  the  two-dimensional  Fourier  trans- 
form yields  the  frequency  wavenumber  power  density  (FWPD)  spectrum. 

Thus,  if  Aj (w)  and  A^(w)  denote  the  discrete  complex  Fourier  spectrum  at 
radial  frequency  w at  antennas  j and  m then  the  cross  spectrum  between  the 
antennas  j and  m at  frequency  w 


In  general,  a smoothed  version  of  cross-sper ';ral  estima  ;es  Sjjjj(w) 
The  estimated  FWPD  spectrum  is  defined  by 


is  used. 
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(2.1.1J 


where  W is  the  total  number  of  antennas,  r^^'s  arc  the  vector  distances  from  some 
arbitrary  origin,  is  the  wavenumber  vector  and  the  *s  are  suitable  weights 
imposed  on  the  j^*^  antenna. 


In  a previous  report  [/,] , different  methods  of  spectral  analysis  to  extract 
frequency,  phase  and  coherence  information  were  described  in  detail.  It  was 
found  that  the  frequency-averaging  technique  provided  a very  efficient  and 
reliable  means  of  cross-spectral  estimation.  Although  it  is  true  that  the 
maximum  likelihood  method  of  spectral  analysis  provides  an  optimal  estimation 
of  the  spectrum,  it  proves  to  be  very  costly  in  terms  of  comjiutcr  time  and 
memory  requirements  and  the  improvement  obtained  thereof  is  not  substantial. 

As  the  number  of  antennas  is  increased,  the  analysis  soon  becomes  impractical. 
Therefore,  only  the  frequency-averaging  technique  of  estimating  the  cross- 
spectra is  used  in  all  of  the  data  involved  here.  Once  the  cross-spectra  have 
been  estimated,  the  spatial  component  of  the  FWPD  spectrum  can  be  estimated  by 
various  methods.  In  the  subsequent  analysis  three  methods  of  estimation  of 
the  angle  of  arrival  have  been  implemented.  The  three  methods  can  be  delineated 
in  terms  of  the  choice  of  the  weights  Wj  and  are:  (1).  The  Fast  Frequency 

Wavenumber  Analysis  or  the  Conventional  Analysis,  which  is  most  useful  for  an 
initial  diagnostic  scan  in  the  frequency-wavenumber  (f-k)  space.  Wavenumber 
sidelobe  characteristics  can  be  varied  by  proper  choice  of  weights  for  the 
different  cross-spectral  terms.  But  resolution  is  poorer  than  methods  2 and 
3.  (2).  The  minimum  redundancy  version  which  is  only  applicable  to  the  data 
from  a strict  minimum  redundancy  configuration  of  the  antennas.  Sidelobcs 
are  better  but  the  technique  suffers  from  grating  lobes  having  amplitudes 
equal  to  tiie  mainlobc  whenever  the  basic  antenna  spacing  exceeds  half  the 
operating  wavelength.  (.3j . An  adaptive  method  using  tlie  maximum  likcliliood 
principle  offers  higli  resolution  in  spatial  characteristics  and  is  applicable 
to  any  arbitrary  antenna  configuration.  Unlike  the  other  two  methods,  this 
method  assigns  adaptive  weights  to  the  different  cross -spectral  terms  depending 
on  the  look-angle  and  the  spatial  characteristics  of  interference  and  noise. 
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Each  of  these  three  approaches  will  be  briefly  outlined  here.  The  conunon 
part  of  the  analysis  and  computation  to  all  three  methods  is  to  obtain 
suitable  cross-spectral  estimates  between  any  given  pair  of  antennas. 


2 . 2 Estimation  of  Cross  Spectra 

The  recorded  data  consists  of  two  sets  of  256  equally  spaced  complex  valued 
time  samples  from  each  antenna  output.  Although  the  entire  256  samples  could 
be  used  for  further  processing,  it  was  decided  to  treat  only  the  middle  128 
samples  in  each  set  for  reasons  of  economy.  Thus,  the  truncated  sets  consisting 
of  128  samples  in  each  subcase  are  immediately  transformed  into  the  frequency 
domain  via  a Fast  Fourier  Transform.  The  resulting  128  complex-valued  discrete 
Fourier  spectra  stretches  across  a basic  Doppler  frequency  range  of  8 Hz.  The 
two  sets  of  time  samples  are  displaced  from  each  other  by  an  amount  g = of 
the  sampling  period.  Therefore,  it  is  possible  to  either  align  or  interlace 
the  resulting  two  sets  of  Fourier  spectra.  If  the  Fourier  spectra  are  aligned 
the  signal -to-ratio  can  be  improved  but  the  basic  Doppler  frequency  range  is 
still  8 Hz.  On  the  other  hand,  if  the  two  setj  are  interlaced  the  frequency 
range  can  be  doubled  to  16  Hz  without  any  improvement  in  the  signal -to-noise 
ratio. 

In  the  case  of  signals  transmitted  by  an  aircraft,  the  basic  range  of  8 Hz 
is  sufficient.  But  ground  bacKscattjr  signals  often  occupy  a band  wider  than 
8 Hz,  in  which  case  interlacing  proves  to  be  very  useful.  In  any  case,  the 
combined  or  interlaced  Fourier  spectra  from  the  starting  point  for  further  data 
processing. 

Three  types  of  output  have  been  designed.  These  consist  of  (1)  power 
spectra,  (2)  coherence  spectra,  and  (3)  maps  of  Doppler  frequency*vs*arrival 
angle  called  DAASM  maps. 


2 . 3 Power  Spectra 


Power  spectra  are  calculated  for  each  antenna  by  taking  the  product  of  eacii 
I'ouricr  spectral  line  with  its  own  complex  conjugate.  Of  course,  the  Fourier 

\ 
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spectra  are  normalized  in  each  channel  so  that  the  total  power  in  each  channel 
equals  unity.  That  is,  if  represents  the  unncrmalized  frequency  line, 
then  the  normalized  Fourier  spectrum  is  simply 


g 


mA 


(2.3.1) 


The  power  spectri  are  smoothed  by  the  frequency-averaging  technique  [2]  ove^ 

7 frequencies  concurrently  with  the  operation  of  evaluating  the  product  of  the 
Fourier  spectra.  This  feature  enables  the  cross-spectra  to  he  estimated  in  an 
identical  way  merely  by  inserting  the  Fourier  spectra  of  the  appropriate  pair  of 
antennas.  The  power  spectra  and  the  cross-spectra  are  different  only  in  a 
matter  of  identification  of  antenna  elements. 


The  power  spectra  thus  obt.iined  are  plotted  on  a standard  digicoder  format 
of  logarithmically  scaled  numbers  ranging  from  1 to  15  which  corresponds  to  a 
linear  scale  from  0 to  100.  TTie  exact  equation  defining  the  scales  is 


DI  = 2 ♦ 7 logjQ  P lin 


1 < P, . <100 

lin 


= 2 X P 


lin 


lin 


< 1 


where  UI  is  the  digicoder  output  integer  and  •’lin  is  the  value  of  the  power 
designed  to  have  a maximum  of  100. 

The  details  of  the  Preface  and  Data  format  are  given  in  [4].  The  power 
spectral  digicoder  printouts/outputs  provide  the  basis  for  selection  of  data 
segments  suitable  for  more  complex  and  time  consuming  analysis.  A representative 
segment  of  data  is  presented  in  Figure  11.  The  data  as  usual  consists  of  sub- 
cases where  each  subcase  has  four  different  height  gates  and  two  carrier 
frequencies.  Therefore,  for  the  six  antenna  configuration  there  will  be  48 
channels  of  power  spectra  for  each  subcase.  A preface  line  separates  each  height 
gate  and  the  subcases  themselves  are  separated  from  each  other  by  6 preface 
lines . 

j 
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2.4  Coherence  Spectra 


Coherence  spectra  between  the  and  n^h  channels  is  defined  by  the 
equation 


coh(f) 


1 , 

Js  (f)  S (f) 
f mm  nn'- 


(f  =,— ) 

2 n 


(2.4.1) 


A- 

where  S^j^,  is  the  smoothed  cross  spectrum  at  frequency  f between  m^*^  and  n^h 
channels  and  S,„j„(f)  would  be  the  self  power  spectrum  of  the  mth  channel  and  so 
on.  It  is  clear  that  the  coherence  is  a number  which  lies  between  0 and  1. 

This  linear  range  is  converted  to  an  arc  hyperbolic  tangent  scale  and  multiplied 
by  a suitable  scale  factor  to  provide  the  digicoder  output  numbers. 


Table  12  shows  the  relationship  between  coherence  and  digicoder  output  num- 
bers. The  coherence  spectra  are  plotted  for  each  available  distance  between  the 
various  pairs  of  antennas  in  order  of  increasing  distance.  Since  the  coherence 
is  only  available  at  specified  distances  the  plots  are  backfilled  at  intermediate 
distances.  It  was  found  the  raw  coherence  thus  plotted  appeared  sometimes  to  be 
a discontinuous  function  of  distance  instead  of  raonotonically  decreasing  fur.ction. 
Therefore,  three  successive  subcases  have  been  averaged  for  each  distance  at  each 
Doppler  frequency  to  provide  a much  smoother  picture  and  thus  the  statistical 
confidence  limits  which  depend  on  any  given  sample  size  are  improved  (see  Fig.  11). 


The  coherence  is  a good  measure  of  the  spatial  extent  of  the  frequency- 
wavenumber  spectruiii,  when  the  number  of  sensors  is  small.  Looking  from  another 
point  of  view,  the  coherence  can  be  interpreted  as  a measure  of  the  relative 
amounts  of  propagating  and  non-propagating  noise.  However,  direct  evaluation 
of  the  FWPD  is  a superior  approach  and  contains  all  the  information. 

The  coherence  spectra  are  primarily  useful  to  determine  the  maximum  meaning- 
ful aperture  in  terms  of  operating  wavelength  for  the  type  of  signals  involved. 

A smooth  ionosphere  (certain  ionospheric  modes)  produces  signals  with  high 
spatial  stationarity . In  such  a case  the  coherence  does  not  decay  rapidly  with 
Increasing  distance.  On  the  other  hand,  a turbulent  ionosphere  displays  high 
coherence  over  relatively  smaller  apertures.  Digicoder  plots  of  the  coherence 
spectra  are  presented  for  various  selected  data  segments  elsewhere  [4] . 
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2.5  Frequency-Wavenumber  Spectra 


DAASM  maps  consist  of  digicoder  plots  of  frequency  against  wavenumber 
which  ii:  turn  is  easily  related  to  the  angle  of  arrival. 


A particularly  efficient  algorithm  to  calculate  the  f-k  spectra  is  possible 
if  the  weights  in  Equation  (2.1.1)  are  chosen  to  correspond  to  a triangular  taper 
function  and  will  be  referred  to  as  Fast  Frequency  Wavenumber  analysis  or  the 
generalized  Bartlett  estimate.  In  such  a case,  the  individual  cross -spectral 
terms  need  not  be  separately  calculated. 


'fhe  separate  Fourier  spectra  from  each  antenna  output  can  be  used  directly 
and  frequency  smoothing  can  be  applied  in  the  final  step.  The  details  of  the 
algorithm  are  as  follows: 

Let  |Aj  (w)  I exp  i(J>j  (w)  represent  the  Fourier  spectral  line  at  frequency  w 
from  the  jth  antenna.  Then  the  cross-spectrum  between  the  jth  and  mth  antennas 
is 


Sjm(w)  = |Aj(w)|  |Aj|j(w)|  exp  i(4»^  - (J)J 
and  the  f-k  spectrum 

N N 

P(w,it)  = s.^(w)  exp 

j=l  m=l 

where  it  is  the  wavenumber  and  where  w = 2Tif. 


(2.5.1) 


(2.5.2) 


Equation  (2.5.2)  can  be  written  alternately  as 

N 

P(w.it)  = |A^(w)|  exp  ij<|>j(w)  - icT^ 

j = l 
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The  subscripts  j and  m can  be  replaced  by  a single  subscript  n and  the 
double  siiinination  can  be  replaced  by  a single  sunmation.  Thus, 


P(w.k)  = 


E 

n=l 


A^(w)  exp  i - k*r^j 


(2.5.3) 


Equation  (2.5.3)  can  be  computed  much  faster  than  equation  (2.5.2). 

P(w,k)  is  calculated  at  as  many  points  as  desired  in  k,  say  with  increments 
Ale.  Each  successive  point  in  Ic  can  be  computed  from  the  previously  computed 
point  by  simple  complex  multiplication  because  exp  i (k  + Ak)  = exp  ik*exp  iAk. 
This  step  is  common  to  all  three  approaches  of  f-k  spectral  calculations  and 
saves  considerable  computer  memory  and  time. 

Since  frequency  smoothing  is  independent  of  the  operations  on  ^ a set  o. 
weights  W(w)  is  calculated  using  P(w,o)  for  each  w and  each  P(w,k)  is  multiplied 
by  the  appropriate  W(w)  to  yield  P(w,k),  the  frequency  smoothed  estimate.  The 
relation  between  the  wavenumber  k and  the  direction  of  arrival  9 for  the  case  of 
a linear  array  is  given  by 


I . |k|  = 

where  k is  now  a scalar,  f the  operating  carrier  frequency  and  c the  velocity 
of  light. 

The  minimum  redundancy  version  can  be  treated  as  a special  case  of  the 
general  equation  for  the  f-k  spectra  given  by  Equation  (2.1.1).  It  is  assumed, 
of  course,  that  a certain  number  of  consecutive  multiples  of  the  basic  antenna 
spacing  are  available  by  suitably  placing  the  different  antennas.  Distances 
beyond  the  largest  consecutive  multiple  are  regarded  to  have  zero  weight,  i.e., 
their  contribution  to  the  f-k  spectra  is  nil.  The  non-zero  weights  decrease 
as  a cosine  squared  function  with  distance.  In  addition,  the  self-spectral  term 
is  assigned  a weight  of  0.5,  so  that  the  expression  for  the  f-k  spectrum  can  be 
written  as 
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L 

P(w,it)  = 1/2  exp  i (k.dp 

Jl=l 

where  Sq  corresponds  to  the  self  term  and  is  the  cross-spectral  term  corres- 
ponding to  the  pair  of  antennas  separated  by  a distance  dj^.L  is  the  largest 
successive  multiple  of  the  basic  spacing  and  is  the  weight  given  by  the 
cosine-squared  function.  It  is  well  known  that  the  cosine-squared  weights  result 
in  good  sidelobe  suppression.  It  can  be  shown  that  for  six  antennas  only  ten 
consecutive  multiples  are  possible  and  for  eleven  antennas  the  corresponding 
number  is  45.  A similar  configuration  for  more  than  eleven  antennas  is  not 
known  at  present.  Each  distance  occurs  only  once  and  hence  ^he  name  'minimum 
redundancy'  for  the  array.  Only  those  cross -spectral  terms  corresponding  to  the 
multiples  of  the  basic  spacing  need  be  calculated.  Thus,  for  eleven  antennas 
45  of  66  possible  terms  are  necessary  resulting  in  almost  a 30%  reduction  in 
calculations.  More  than  that,  in  a sense,  the  weights  are  almost  optimal  for 
the  configuration  and  there  is  excellent  sidelobe  suppression.  The  carrier 
frequency  can  be  chosen  such  that  the  basic  spacing  is  a half-wave  length  or 
smaller,  thus  avoiding  grating  lobes. 

The  maximum  likelihood  FWPD  spectrum  is  derived  from  sets  of  filters  which 
are  data-adaptive  and  optimal.  The  filters  are  optimal  in  the  sense  that  a unit 
plane  wave  propagating  across  the  array  of  sensors  arriving  from  a particular 
direction  is  passed  undistorted  while  signals  and  noise  from  other  directions  are 
minimized  in  a mean  square  error  sense.  When  the  noise  distribution  is  assumed 
to  be  multi -dimensional  Gaussian,  the  results  are  identical  with  the  maximum 
likelihood  estimate  of  the  signal  calculated  from  the  likelihood  functions  of 
signal  plus  noise. 

This  method  of  f-k  r.pectral  estimation  employs  ar  adaptive  scheme  to  gene- 
rate the  set  of  weights  in  the  filter  given  by  Equation  (2.1.1).  Instead  of 
using  a fixed  set  of  weights  for  all  k the  adaptive  schem.e  allows  for  the  set 
of  weights  to  be  varied  depending  on  the  particular  k as  well  as  signal,  inter- 
ference and  noise  characteristics.  In  other  words,  the  set  of  weights  is  dictated 
by  the  data  sample  available,  subject  to  the  constraints  that  a unit  plane  wave 
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signal  travelling  with  a given  wavenumber  k = Kq  (arriving  from  the  direction 
towards  which  the  antenna  beam  is  pointed),  is  passed  undisturbed  through  the 
spatial  filter  but  everything  else  for  which  k ^ Kq  (arriving  from  other  direc- 
tions) is  rejected  in  an  optimal  manner.  Thus,  the  true  spectrum  is  reproduced 
with  high  resolution  and  not  smeared,  which  is  the  result  of  convolution  between 
the  true  spectrum  and  the  particular  spectral  window  chosen.  This  choice  of 
constraints  results  in  a convenient  mathematical  formulation.  It  is  only  neces- 
sary to  calculate  the  inverse  of  the  matrix  of  cross -spectral  terms  for  each 
Doppler  frequency  of  interest,  in  order  to  furnish  the  optimal  set  of  weights 
automatically  adapted  to  each  k for  a given  data  sample. 
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3.  MATHEMATICAL  DERIVATIONS  AND  PROPERTIES  OF  THE  MLM  ESTIMATOR 


3.1  Minimum-Variance  Unbiased  Estimate 


Consider  the  set  of  measurements 

z 

where  z is  the  measurement  vector, 
additive  random  noise  vector  and  H 
set  that  minimizes  the  mean  square 

A 

estimated  vector  x then 

where  the  superscript  H is  the  conjugate  transpose.  If  A is  so  chosen  that  v is 
rejected  by  the  filter  but  Hx  is  passed  undistorted  then 

Hx  = Ha”z  = HA^Slx 

H H 

or  1 = HA  and  by  symmetry  AH  = 1 which  is  the  constraint  expressed  mathematically. 

A 

Now  the  variance  of  x can  be  written  as 

= a“rA 

where  R is  the  covariance  matrix  of  the  noise  process. 

If  is  minimized  under  the  constraint 


-*■  -+• 

= Hx  + V 

X is  a constant  signal  vector,  v is  the 
is  the  signal  model  matrix.  If  A is  the  filter 
error  between  the  signal  vector  x and  the 


h“a  = 1 


Then  a quadratic  cost  function  J can  be  written  as 

J = A^RA  + X(Ha“-1) 


J is  minimized  when 


0 
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I 


or 


a“r  + xh“  = 0 

A»  = -XH«R-' 
A^  = 


-A  = + 


a“h  _ 

h“r"^h  h”r‘^h 


= +- 


Ii"r-1 

A = 


h’^r'^h 


= A^z  = (h“r"^H)"^  h“r’^z 


3.2  Properties  of  the  Minimum  Variance  Estimator 


From  Equation  (3.1.1) 


(h»r-1h)  h«r-‘? 


' < '^1!  "k  "kJ  'k 


vviiere  k = dimension  of  the.  measurement  vector  premultiplying  the  equation 


"k  = \ 


by  (h“  R^j  H^)  h|J  R'J  , one  obtains 


-1 


-1 


< \i  »k)  t ^ - '"k  V "k  \l  »k* 


-1 


<"k1v 


(3.1.1) 
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3.3  Maximum  Likelihood  Estimate 

The  previous  derivation  minimizes  the  mean-square  error  subject  to  the 
constraint  that  a unit  plane  wave  be  passed  undistorted  through  the  filter.  But 
there  is  no  probabilistic  basis  for  the  characteristics  of  the  noise  field.  In 
other  words,  the  least  squares  estimate  does  not  assign  probability  density 
functions  to  the  signal  or  measurement  vectors.  The  derivations  are  obtained 
via  deterministic  arguments  only. 

It  is  possible  to  establish  probabilistic  arguments  if  the  random  process 
involved  is  treated  as  multi -dimensional  Gaussian.  The  corresponding  likelihood 
function  can  be  maximized  and  thus  one  arrives  at  the  maximum  likelihood  prin- 
ciple. The  estimated  signal  vector  maximizes  the  probability  of  the  measurement 
vector  that  actually  occurs,  taking  into  account  the  properties  of  the  additive 
random  noise  vector. 

Consider 

->• 

z = Hx  + V 
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where  z is  the  measurement  vector,  x is  the  true  signal  vector  and  v is  the 
random  error  vector  and  H is  an  appropriate  dimensional  matrix  representing  the 
transformation. 

The  likelihood  function 


L(z.x)  = P(z'x) 


P(z>x)  ^ P(X,VJ  3z 
P(x)  P(x)  3^ 


But  — = 1.  Therefore, 
3x 


P(x,y)  3z  _ P(x,v)  _ P(x)  P(v) 
3x  P(x)  P(x) 


= P(v)  ^ Py(l-Hx) 


3L 


P^(z-Hx) 

3x 


0 


The  above  shows  that  the  likelihood  function  L(z,x)  which  is  defined  as 
the  probability  density  function  of  z conditioned  on  x is  just  the  density 
function  of  v centered  around  Hx. 

If  V is  taken  as  a zero-mean-Gaussian  distributed  vector  with  covariance 
matrix  R,  then 


P(zlx)  Tjj  exp  [-^(z-Hx)^  R’^(z-Hx) 

Maximizing  P(?|jt)  is  accomplished  by  minimizing  the  exponent  in  brackets 
with  respect  to  x;  that  is  J = (1-Hx)**  R~^(?-Hx)  where  the  superscript  H denotes 
complex  conjugate  transpose. 


-h“r"^z  + H*^R"^Hx  = 0 

X = (h“r"^)"^  h“r‘^z  (3.3.1) 
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and  is  identical  to  the  weighted  least  squares  estimate.  The  equivalence  of 
the  maximum  likelihood  estimate  for  a multi-dimensional  Gaussian  distribution 
of  the  random  noise  vector  and  the  weighted  least  square  estimate  under  the 
constraint  that  a unit  plane  wave  be  passed  undistorted  through  the  filter  has 
thus  been  established. 


3.4  Estimation  of  MLM  FWPD  Spectra 

The  previous  discussion  showed;  how  in  general,  from  a noisy  measurement 
vector,  an  estimate  of  the  desired  signal  vector  can  be  obtained  such  that  the 
sura  of  the  squares  of  the  deviations  is  a minimum  under  the  "distortion  free" 
constraint.  The  details  of  application  of  the  least  squares  or  maximum  likeli- 
hood principle  to  the  estimation  of  the  FKPD  spectrum  are  set  forth  in  the 
following. 


A unit  amplitude  signal  model  is  represented  by 


i2Tric»r  j! 

i2frk»r2 

» 

i2iric»r 

n 


where  n is  the  number  of  sensors  in  the  array  r^j  is  the  vector  distance  to  the 
n^b  sensor  from  an  arbitrary  origin,  and  Ic  is  the  wavenumber  vector  correspond- 
ing to  the  direction  of  arrival  of  the  signal.  Monochromatic  operation  is  assumed. 


If 
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} 


is  a filter  set,  then  the  mean-square  error  for  this  filter  sot  is 

e » (1-f“e)  (1-e“f)  + f“sf 

e X i_e“f  - f“e  + f“ee“f  + f“sF  (3.4.1) 

where  S is  the  noise  cross  spectral  matrix  and  the  sxoperscript  H denotes  the 
conjugate  transpose. 

The  mean-square-error  (MSB)  should  be  minimized  subject  to  the  constraint 

that 

f“ee“f  1 

since  this  represents  the  filter  power  output  for  a unit  amplitude  plane  wave 
signal  model. 

This  constraint  can  be  incorporated  if  the  expression 

e » (1-e“f)  (1-f“e)  + f“sf  + X(f”ee“f-1)  (3.4.2) 

is  minimized,  where  X is  the  Lagrangian  multiplier.  X c?ji  be  determined  by 
taking  partial  derivatives  with  respect  to  the  real  and  imaginary  components  of 
F. 

Consider  the  m^^  component  of  F 

^ - E„a-E”F)  . (SF)„ 

m 

+ (f“s)„  + x|e*(f“e)  + E (e“f) 
m I m ***  J 

and 
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But  for  an  extremum 


j Se  _ 

5ReF“  5lmF“ 

m m 


. e = -2E  (1-e“f)  + 2(SF)  + 2XE  (f“f) 

m mm 


SF  = E(1-E^F  - XE*^F) 


F = [1-(A+1)  e“f]S‘^E 
F^EE^F  = [f“e-(X+1)  (F^^Ee’^F)]  e”s-1E 


Since  F^EE^F 


= 1 (by  the  constraint) . 


-f^^=  f“e-(X+i) 
E S E 


^ 1 rHr  1 
A+1  = F E + — r= = 

e“s"  E 


SF  = eTi  - (f“e  + eM 

L E S~^E  J 


e"s’^e 


_ F^EE^F  _ 1 

r or  - u -I  ” li  i 

E S"  E E S"  E 


A 


Therefore,  the  expression  for  the  MSE  becomes 


e = 1 


-E*V| 


2 » 1 


and  is  minimized  obviously  when  E^F  = 1 = a condition  consistent  with 

f‘Wf  =1. 

Since 


SF  = 


ee”f 

E^S'^E 


from  Eq.  (3.4.5)  the  filter  set  which  minimizes  the  MSE  is 


F = S‘^E/e“s-1E 


(3.4,6) 


and  the  minimized  noise  power 


f“sf  = 


e“s-ie 


(3.4.7) 
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Also  e F where  0 is  an  arbitrary  phase  factor  accomplishes  the  same  task. 


The  filter  F also  maximizes  the  signal  to  noise  ratio  G = 


F^EE^F 

f“sf 


without  being  subjected  to  the  constraint  that  F^EE^F  = 1.  Similarly  Ce^®F 

10 

would  maximize  the  signal  to  noise  ratio  where  Ce  is  an  arbitrary  complex  valued 
scalar. 


3 . 5 A Relevant  Theorem  and  Additional  Properties 


The  foregoing  can  be  neatly  summarized  as  a theorem  involving  Hermitian 


forms . 


G = 


F^SF 


(5.5.1) 
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If  a quantity  G is  expressed  as  a ratio  of  two  Hermitian  forms  where  S is 
nonsingular  and  positive -definite,  the  largest  eigenvalue  of  the  regular 
pencil  of  matrices  [U-XS]  is  the  maximum  obtainable  G when  F is  the  eigen- 
vector satisfying  the  homoegenous  equation 

UF  = Xj^^F 

If  U is  expressible  as  U = EE^  then  the  largest  and  non-zero  eigenvalue  of  the 
regular  pencil  [U-XS]  is 


(3.5.2) 


and  the  eigenvector  corresponding  to  X,,  is 

M 


(3.5.3) 


The  above  equations  solve  the  problem  of  optimization  for  spatial  process- 
ing of  array  data. 

The  Hermitian  form  E^*SE  can  be  reduced  by  a unitary  transformation  E = Tx 
to  the  canonical  form 


N N 

^A.XjXi  =Vxj|xj 
i=l  i=l 


(3.5.4) 


where  the  -Vs  are  the  eigenvalues  of  S and  the  x^’s  are  the  components  of  E 
along  the  eigenvectors  of  S. 

Therefore  the  eig«ivalues  of  S"^  would  be 


%’X,”  • 


and 


P 


L 


E^S'^E 


where  is  the  maximum-likelihood  filter  output  power. 


(3.5.5) 
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If  the  probe  vector  E lies  along  the  eigenvector  then  = Xj  and  reflects, 
the  correct  value  of  the  spectrum. 

If  is  the  conventional  filter  power  output,  the  relative  resolution  can 
be  found  from  the  ratio 

L 1 3 X. 


= E Z 
i j 


2 


Therefore, 


since 


(sUiP)’ 


X. 

* 1 z ^ 

3 


(is  positive  real) 


1 


• 1 


(3.5.6) 


Thus,  the  MLM  spectrum  has  lower  output  power  off  the  principal  axes  of 
the  spectral  matrix  S and  yet  Pq  = P^  along  the  principal  axes.  The  MLM  spectrum 
exhibits  sharper  peaks  and  higher  resolution.  Although  this  result  is  not 
strictly  applicable  when  the  eigenvectors  of  S lie  outside  the  space  spanned  by 
the  subset  of  pi"»be  vectors  corresponding  to  wavenumber  vectors  if,  in  practice 
the  result  quite  robust. 


3.6  Random  Noise  and  TVo  Plane  Waves 

The  behavior  of  the  maximum  likelihood  FWPD  spectrum  can  be  analyzed  when 
the  noise  field  consists  of  a small  number  of  unit  plane  waves  and  uncorrelated 
noise  propagating  across  an  array  of  sensors. 
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The  analysis  is  based  on  the  assumption  that  the  cross  spectral  matrix  of 
such  a noise  field  can  be  represented  as 

N 

S = PI  +7^ V v”  (3.6.1) 

n n 

n=l 

where  I is  the  unit  matrix  and  P the  variance  of  the  noise.  V is  the  n^^  unit 

n n 

plane  wave  dyadic  representation.  This  representation  is  a matrix  of  unit  rani;. 

It  does  not  have  much  physical  significance  when  a large  number  of  plane  waves 
are  present. 

Consider  the  case  of  two  complex  sinusoids  and  white  noise.  Then 

s = PI  + uu“  + w“ 


where 


i2Trk 


U = a X col(e 


A « 

u 1 


i2iTk 


u 


...) 


and 


i2iric  »x,  i2TTk  ‘x- 
V = b X col (e  ^ , e ^ . . . .) 

k^  and  k^  representing  the  wavenumbers  for  the  two  plane  waves. 

Therefore,  the  MIM  spectrum  can  be  shown  to  be 


p 

1 

P+|U| 

V« 

^ i/v 

p+|v|^ 

b 1 

E S"  E 

|E|^ 

e"u 

E»V 

V*^E 

P+|u|* 

A 

V*^E 

v»u 

p+|v 
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1 


= 0 if  Ej<au-*-6V  where  a and  3 are  complex -valued  scalars. 


from  which  lim 

p-»o  e"s’*E 


In  particular  if  E = IJ/|U 


and  if  E = V/|v| 


lim 


p-*o 


E*^S‘^E 


|U 


2 


lim 

p-KJ 


1 

E^^S'^E 


V| 


2 


Therefore,  the  spectrum  exhib  ti»  i type  of  pseudolinearity  in  the  vicinity 
of  the  directions  of  U and  V provided  P is  small.  The  conventional  estimate 
E^SE  yields 


lim  e”sE  = |U|^  + I-— i ^ 

P-o  |U|^ 


if 


and 


if 


Iv'l 


The  MLM  spectrum  gives  the  true  value  of  U and  V when  the  steering  vector 
points  in  either  of  these  dire(,tions  whereas  the  conventional  estimate  is  corrupted 
by  the  influ'?nce  of  U while  pointed  towards  V and  vice  versa.  At  all  other  points 
E,  not  a linear  combination  of  U and  V,  the  spectrum  tends  to  zero  when  P is 
small . 
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In  the  case  of  a single  con^lex  sinusoid  if 


then 


lim 

p-^ 


e”s'^e 


and 


lim  e“sE  = |v|* 

p-K) 


whereas  if  E aV 


lim  -T.— -H--  » 0 
p*o  E S"  E 

but 

lim  e“se  = lu“v|^ 

p-K) 

reflecting  the  window  effect. 


3.7  Computational  Procedure 

The  computational  procedure  is  set  forth  in  the  following  steps. 

Hie  normalized  input  Fourier  spectra  from  the  various  antennas  is  converted 
to  cross-spectral  estimates  Sj^^^Cw)  by  anyone  of  several  spectral  estimation  pro- 
cedures, in  this  case  by  frequency  averaging  and  smoothing  with  seven  frequencies. 
Subscripts  j and  m refer  to  the  antennas  in  question.  Next  an  NxN  matrix  of 

A 

cross -spectral  estimates  [S(w)]  is  organized  at  each  w,  where  N is  the  total 
number  of  antennas.  Not  all  NxN  terms  need  be  evaluated  since  Sjjj(w)  = Sjjj(w), 
where  the  asterisk  denotes  complex  conjugation.  The  sum  of  the  diagonal  elements 
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in  the  matrix  is  stored  as  a measure  of  total  power  at  any  given  Doppler 
frequency  and  used  for  renormalization  purposes  ultimately.  The  matrix  itself 

is  normalized  by  replacing  each  element  Sjj,(w)  by  Sj„(w)^ySjj  (w)Sj^(w) , a step 

which  ensures  sensor  equalization.  The  normalized  cross-spectral  matrix  is 
inverted  for  each  Doppler  frequency  one  at  a time.  The  elements  of  the  inverted 
matrix  are  used  in  the  equation  for  f-k  spectra 


P(w.ic) 


r N N 


m=l  n*l 


exp  ik*(r.-r  ) 
j In 


1 


It  is  seen  that  the  weights  w^^  and  the  straightforward  cross -spectral  elements 
together  are  replaced  by  the  corresponding  inverted  matrix  elements  and  the 
reciprocal  is  calculated  to  yield  the  high-resolution  f-k  estimate.  The  incre- 
ments over  k are  calculated  in  a chain  sequence  from  the  previous  values  in  the 
same  manner  as  indicated  before.  The  final  step  is  to  restore  relative  power 
levels  at  the  different  Doppler  frequencies.  First,  the  values  of  P(w,k)  are 
summed  over  all  k and  the  normalized  values  of  P(w,k)  are  multiplied  by  the  sum 
of  the  diagonal  elements  PjsC'*)  original  cross -spectral  matrix  which  have 

been  stored  as  indicated  before.  Thus,  the  renormalized  P(w,k)  are  obtained 
and  with  an  appropriate  scale  factor  and  conversion  to  a digicoder  format  are 
recorded  on  tape  to  produce  the  digicoder  printouts. 
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4.  RESULTS 


4.1  Simulated  Models  of  Signals  and  Random  Noise 

This  section  discusses  the  results  of  the  application  of  the  maximum 
likelihood  method  to  simulated  models  of  signals  and  white  noise  generated  in 
complex  valued  digital  form  on  the  computer.  The  signal  part  of  the  test 
function  is  modeled  as  the  sum  of  one  or  more  complex  sinusoids  with  appropriate 
progressive  phase  shifts  along  the  linear  array.  Thus,  if  two  sources  on  two 
different  Doppler  angular  frequenc  .es  Wj  and  Wj,  and  located  in  directions  9i 
and  02  from  the  axis  of  the  array  have  to  be  simulated,  the  expression 

Sjj(t)  = Aj  exp  j(wjt  ♦ ^ cos  6,) 

+ Aj  exp  j (W2t  + djj  cos  02 ) 

is  used.  Here  Aj  and  Aj  are  the  respective  amplitudes  of  the  two  sinusoids,  X 
is  the  operating  wavelength  of  the  r*f  carrier,  and  the  subscript  n is  the 
antenna  index,  d^  being  the  distance  of  the  n^*'  antenna  from  the  origin.  The 
value  of  s„(t)  is  calculated  at  discrete  equally  spaced  time  points  t^,  t2»  etc., 
to  obtain  the  sampled  data  points. 

Additive  white  noise  is  introduced  by  appropriately  scaled  random  numbers 
from  a random  function  generator.  The  scaling  determines  the  amount  of  noise 
power  added.  The  sampled  data  points  along  with  the  noise  is  processed  as 
described  earlier  to  produce  DAASM  maps. 

4.2  Parameter  Variations  and  Digicoder  Displays 

The  performance  of  the  MIAi  was  studied  with  test  functions  designed  to  vary 
the  following  parameters. 

1)  Several  noise  levels  (Figs.  1-10) 

2)  Different  number  of  Doppler  frequencies  being  smoothed  (Fig.  1) 

3)  Numbe-  of  sources:  a)  coherent  (same  Doppler)  (Figs.  8 8 10) 

b)  incoherent  (separated  Dopplers  (Figs.  2,  3,  4,  5, 

6,  7 § 9) 
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4)  Proximity  of  sources  a)  wavenumber  (Figs.  4,  5) 

b)  Doppler  frequency  (Figs,  2,  3,  9) 

5)  Number  of  antennas  either  6 or  12  § configuration 

a)  uniform  (Figs.  1,  2,  5,  4,  5,  6,  7) 

b)  minimum  redundancy  (Figs.  6,  7,  9,  10) 

Salient  resul.j  pertaining  to  the  above  variations  are  presented  here  in 
the  form  of  digicoder  plots  of  Doppler  frequency  against  wavenumber  (angle  of 
arrival).  Tables  1 - 10,  provides  the  details  of  the  parameters  in  each  plot. 

The  digicoder  plot  is  designed  so  that  there  are  128  Doppler  frequency 
lines  and  128  wavenumber  lines.  Since  the  width  of  the  digicoder  printout  can 
accommodate  256  lines  two  different  cases  can  be  presented  on  the  same  plot, 
one  in  the  left  half  and  the  other  in  the  right  half. 

The  characters  represent  numbers  from  0 to  15  in  increasing  order  of 
intensity  and  as  described  before  logarithmically  scaled  from  the  corresponding 
linear  values  of  the  spectra.  In  effect,  the  digicoder  plot  is  also  an  intensity 
plot.  Sixteen  preface  lines  precede  the  actual  plot  and  are  utilized  to  code 
all  relevant  information  about  the  parameters.  Also,  this  space  is  divided  into 

16  compartments  across  the  full  font  for  viewing  ease.  A similar  compartKifr.ta- 

( 

lization  is  effected  along  the  length  of  the  page  by  interrupting  the  plot  with 
a preface  line  every  16  lines. 

In  the  plots  presented  here  the  wavenumber  axis  is  across  the  width  of  the 
page  and  Doppler  frequency  along  the  length.  The  wavenumber  axis  corresponds 
to  an  azimuthal  direction  of  arrival  between  -90®  and  +90®  with  0®  corresponding 
to  boresight.  Actually  any  sector  of  the  azimuthal  coverage  can  be  expanded 
to  occupy  the  full  128  lines.  Similarly,  the  middle  of  the  Doppler  axis  ' 
represents  0 Hz  becoming  increasingly  negative  towards  the  top  and  increasingly 
positive  towards  the  bottom.  Once  again,  any  portion  of  the  Doppler  spectrum 
with  any  chosen  center  frequency  can  be  accommodated. 

It  should  be  noted  that  the  wavenumber  axis  represents  the  sine  of  the 
angle  away  from  boresight.  Therefore,  the  angles  corresponding  to  the  horizon 
are  compressed  although  approximately  linear  in  the  region  of  boresight. 
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5.  CONCLUSIONS 


5.1  Incoherent  Sources 

Computations  have  been  made  with  the  noise  levels  corresponding  to  power 
S/N  ratios  of  20  dB,  10  dB,  0 dB,  and  -5  dB.  In  the  case  of  incoherent  sources 
assuming  the  large  aperture,  the  estimate  is  very  accurate  at  low  noise  levels 
and  the  accuracy  is  maintained  even  at  angles  close  to  the  horizon  as  seen  in 
Figs.  6,  7 § 9d,  e § f.  As  long  as  the  sources  do  not  overlap  in  Doppler  for  a 
specified  smoothing  they  can  be  resolved  to  within  1®.  As  the  noise  level 
is  increased,  this  resolution  cannot  be  maintained  and  slight  shifts  in  angular 
position  are  observed.  The  valleys  between  the  peaks  do  not  quite  bottom  out 
to  the  low  value  reached  at  angles  away  from  the  source  locations.  Also,  as 
the  noise  level  is  increased  apart  from  the  noise  appearing  randomly  at  all 
Dopplers  and  all  angles,  the  angular  side  lobes  on  the  signal  Doppler  frequency 
form  a ridge.  The  ridge  is  particularly  obvious  when  the  line  source  happens 
to  be  off  the  lines  on  the  MLM  spectra  actually  computed  and  is  attributable 
to  the  procedure  for  renormalization  of  power  across  the  Doppler  frequencies 
as  seen  in  Fig.  4.  The  ridge  effect  vanishes  when  the  source  is  distributee 
over  wavenumber,  since  in  that  instance,  no  single  line  contributes  inordinately 
to  the  total  power  in  the  Doppler  line  and  hence,  the  actually  computed 
adjacent  lines  are  as  effective. 

There  is  a definite  relationship  between  the  number  of  frequencies  that 
are  used  for  smoothing  and  the  number  of  antennas  (sensors)  employed.  In 
general,  no  problems  are  encountered  in  the  way  of  singular  spectral  matrices 
if  the  number  of  frequencies  exceeds  twice  the  number  of  sensors.  But  for 
longer  arrays  this  procedure  is  impractical.  An  alternate  solution  to  avoiding 
singular  spectral  matrices  would  be  to  add  a small  correction  term  A to  the 
diagonal  terms  of  the  normalized  cross-spectral  matrix. 

The  above  implies  that  resolution  in  Doppler  has  to  be  compromised  to 
achieve  improved  angular  resolution.  As  one  increases  the  number  of  spatial 
points  at  which  the  incident  noise  field  is  measured,  one  can  afford  to  ad:^  a 
little  noise  to  the  diagonal  terms.  For  this  reason,  a 15  point  smoothed 
frequency  averaging  is  employed  thus  allowing  adequate  Doppler  resolution 
(Ref.  Figs.  3-10).  In  fact,  it  was  observed  that  5 point  frequency  smoothing 
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resulted  in  ill  conditioned  matrices  for  the  6-antenna  mode  while  7-frequency 
smoothing  was  marginal  (Ref.  Fig.  2).  The  successful  inversion  of  7-frequency 
smoothed  cross -spectral  matrices  can  be  attributed  to  the  noise  added  in  the 
original  time  series.  In  comparison,  the  twelve  antenna  mode  cross -spectral 
matrices  could  only  be  inverted  for  11  or  more  points  smoothing. 

A similar  argument  can  be  advanced  as  to  the  ideal  location  of  the  antennas. 
For  instance,  a small  number  of  antennas  can  be  placed  in  a minimum  redundancy 
configuration  provided  more  time  samples  are  available  to  maintain  a decent 
signal -to-noise  ratio.  A nv *e  reliable  estimate  can  then  be  obtained.  In  a 
certain  restricted  sense,  the  spatial  and  temporal  points  are  interchangeable. 
Conversely,  redundancy  in  'ivailable  distances  between  pairs  of  antennas  in  a 
uniformly  spaced  array  provides  a better  estimate  in  the  presence  of  noise.  If 
angular  resolution  is  the  primary  objective,  it  is  advisable  to  utilize  the 
maximum  possible  aperture  consistent  with  requirements  of  spatial  coherence. 

For  a given  amount  of  frequency  averaging,  this  would  entail  a sacrifice  in 
the  stability  of  the  estimates. 

It  is  observed  that  there  is  coupling  between  the  sources  even  though  they 
are  on  different  Doppler  frequencies  beyond  the  bandpass  edges  of  the  smoothing 
filter.  The  problem  seems  to  be  accentuated  at  low  noise  levels  and  small 
apertures  (Ref.  Fig.  5).  If  the  spillover  pips  are  above  the  background  noise 
levels  the  method  exhibits  ambiguity  as  to  the  location  of  the  sources  cor- 
responding to  a particular  Doppler  frequency.  When  the  background  noise  level 
obscures  the  weakly  coupled  source,  this  phenomenon  vanishes.  The  explanation 
for  the  stronger  coupling  in  the  case  of  smaller  apertures  can  be  traced  to 
the  wider  beamwidth,  and  hence  decreased  resolution. 

However,  this  coupling  phenomenon  does  not  seem  to  be  of  much  concern  in 
the  case  of  actual  measured  data  of  the  DAASM  experiment  because  the  noise 
levels  are  such  as  to  submerge  the  level  of  the  coupled  signal. 

5.2  Coherent  Sources 

Thus  far,  the  results  discussed  pertain  to  the  case  of  incoherent  sources 
and  tnn  technique  suffer*:  from  no  serious  difficulties.  When  the  input  data 
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results  from  coherent  discrete  sources,  the  sidelobes  in  wavenumber  space  reach 
disconcerting  levels  (Ref.  Figs.  8 § 10).  In  general,  it  is  seen  that  two 
sources  of  equal  amplitudes  result  in  an  apparent  sidelobe  level  of  approxi- 
mately -6  dB  or  higher.  This  imposes  quite  a degree  of  uncertainty  in  being 
able  to  distinguish  between  strong  and  weak  sources.  The  sidelobes  could  be 
mistaken  for  weak  sources  or  vice  versa.  Nevertheless,  under  similar 
circumstances  the  conventional  FWPD  spectrum  behaves  even  worse,  true  to  the 
theoretical  predictions  of  the  last  Chapter.  The  ability  to  discriminate 
between  the  two  coherent  sources  also  deteriorates  as  the  noise  increases. 

The  situation  is  improved,  but  not  very  much  when  the  number  of  antennas  is 
increased  from  6 to  12  (compare  Fig.  8 with  Fig.  10). 

The  angular  separation  of  the  sources  beyond  the  natural  beamwidth  for 
any  given  aperture  does  not  seem  to  provide  any  readily  predictable  pattern  of 
sidelobe  location  or  level.  The  source  locations  are  still  discvrnible  with 
a priori  knowledge  but  in  the  case  of  actual  data,  for  instance  coherent 
groundscatter  returns,  it  is  extremely  difficult  if  not  impossible  to  locate 
the  authentic  reflection  points. 

It  is  evident  that  the  case  of  multiple  coherent  sources  has  to  be  solved 
by  either  imposing  additional  constraints  or  by  decoupling  the  sources.  Im- 
posing additional  constraints  would  require  a priori  knowledge  of  the  location 
of  the  sources  which  is  impractical.  In  a real  environment,  the  location  of 
interferers  may  change  or  the  ionospheric  conditions  may  result  in  multipath 
propagation  leading  to  apparent  multiple-source-type  strongly  coherent  signals. 

A method  of  decoupling  the  sources  utilizing  adaptive  feedback  techniques 
can.  be  employed.  An  initial  run  would  yield  the  approximate  location  of  true 
as  well  as  virtual  sources  based  on  a decision  oriented  analysis.  A decoupling 
matrix  can  then  be  constructed  which  in  turn  would  be  incorporated  into  the 
input  cross-spectral  matrix  for  each  wavenumber  of  interest.  The  output  would 
then  conform  to  just  one  source  location  as  though  all  the  other  sources  were 
not  present,  which  should  result  in  very  low  sidelobe  level. 

The  MLM  FWPD  spectrum  has  been  previously  used  in  seismological  applications. 
But  the  application  of  the  technique  to  the  analysis  of  ionospheric  motion  and 
irregularities  is  thought  to  be  new.  It  is  a powerful  technique  with  superior 
resolution,  is  data-adaptive  and  computationally  not  much  more  complex  than 
conventional  methods. 
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TABLE  1. 


Legend  for  Fig.  1 (Effect  of  Smoothing) 


Figs. 
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TABLE  2. 


Legend  for  Fig.  2 (Doppler  Proximity) 


Figs. 
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Antenna  Array;  U = uniform 
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Legend  for  Fig.  3 (Doppler  Proximity) 


Figs. 

Antenna  Array;  U = uniform 

Number  of  Antennas 

Antenna  Aperture 
(Wavelengths) 

Doppler  Lines  (No.) 

Noise  Level  as  % of 
Signal  on  Line  SP 

Doppler  Smoothing,  No.  of 
Points 

Spatial  Angle  of  Sources 
from  Boresight. 
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Legend  for  Fig.  4 (Angular  Proximity,  Incoherent  Sources) 


Figs. 

a 

b 

c 

d 
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Antenna  Array; 
U = uniform 

U 
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U 

U 

U 
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<3 
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<3 

Doppler  Lines  (No.) 

24,59 
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Noise  Level  as  Is  of 
Signal  on  Line  59 

10 

100 

300 

10 

100 

300 
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Legend  for  Fig.  5 (Angular  Proximity,  Incoherent  Sources) 


Figs. 

Antenna  Array; 

U = uniform 

Number  of  Antennas 

Antenna  Aperture 
(IVavelengths) 

Doppler  Lines  (No.) 

Noise  Level  as  % of 
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of  Points 

Spatial  Angle  of  Sources 
from  Boresight. 
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Legend  for  Fig.  6 (Antenna  Configurations) 


Figs. 

Antenna  Array;  U = uniform, 
MR  = minimum  redundancy 

Number  of  Antennas 

Antenna  Aperture 
(Wavelengths) 

Doppler  Lines  (No.) 

Noise  Level  as  % of 
Signal  on  Line  59 

Doppler  Smoothing,  No.  of 
Points 

Spatial  Angle  of  Sources 
from  Boresight. 
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Legend  for  Fig.  7 (Antenna  Configurations) 


Figs. 

Antenna  Array;  U = uniform, 
MR  = minimum  redundancy 

fAunber  of  Antennas 

Antenna  Aperture 
(Wavelengths) 

Doppler  Lines  (No.) 

Noise  Level  as  % of 
Signal  on  Line  59 

Doppler  Smoothing,  No.  of 
Points 

Spatial  Angle  of  Sources 
from  Bores ight. 
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Legend  for  Fi<j.  8 (Coheren':  Sources) 


Figs. 

Antenna  Array;  U = uniform, 
= minimum  redundancy 

Number  of  Antennas 

Antenna  Aperture 
(Wavelengths) 

Doppler  Lines  (No.) 

Noise  Level  as  % of 
Signal  on  Line  59 

Spatial  Angle  of  Sources 
from  Boresight. 
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04 

o 
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TABLE  9. 


Legend  for  Fig.  9 (Antenna  Configuration) 


Figs. 

Antenna  Array;  U = uniform 
MR  - minimum  redundancy 

Number  of  Antennas 

Antenna  Aperture 
(Wavelengths) 

Doppler  Lines  (No.) 

Noise  Level  as  % of 
Signal  on  Line  59 

Spatial  Angle  of  Sources 
from  Boresight. 


a 

b 

U 

U 

12 

12 

<3 

<3 

44,59 

44,59 

10 

100 

0",-30*  0®,-30® 


c 

d 

U 

MR 

12 

11 

<3 

>20 

44,59 

44,59 

300 

10 

0“,-30*  0°,-30“ 


e 

f 

MR 

MR 

11 

11 

>20 

>20 

44,59 

44,59 

100 

300 

0*,-30®  0“,-30® 
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TABLE  10. 


Legend  for  Fig.  10  (Coherent  Sources) 


Figs. 


Antenna  Array;  U » uniform 
MR  = minimum  redundancy 

Number  of  Antennas 

Antenna  Aperture 
(Wavelengths) 

Doppler  Lines  (No.) 

Noise  Level  as  % of 
Signal  on  Line  59 

Spatial  Angle  of  Sources 
from  Boresight. 


a 

b 

U 

U 

12 

12 

<3 

<3 

59 

59 

10 

100 

0",-30*  0*,-30* 


c 

d 

U 

MR 

12 

11 

<3 

>20 

59 

59 

300 

10 

0*,-30‘'  0®,-30“ 


e 

f 

MR 

MR 

11 

11 

>20 

>20 

59 

59 

100 

300 

0*,-30* 

0",-30 
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TABLE  11. 


For  Data  Sample:  Flight  4-196,  15  July  1974,  UT  2100 

See:  DAASM  Map,  p.  98  - ref.  report* 

Coherence  Map,  p.  99  - ref.  report* 


Antenna  Array  - Miniroian  Redundancy,  6 Antennas 

Antenna  Aperture  — 720  Meters 

Doppler  Smoothing,  No.  of  Points  --  15 

Freq  1 (Ffiiz)  --  8.22  Figs,  a,  b,  c 

Freq  2 (MHz)  — 6.62  Figs,  d,  e,  f 

Range  (KM)  — 960 

Azimuth  (deg)  — 349T 

Center  Freq  (Hz)  — +6 


Power  Spectra:  Figures  a 5 d 
Coherence  Spectra:  Figures  b 8 e 
DAASM  Maps,  F-K  Spectra:  Figures  c 5 f 


Sales,  G.S.,  et  al  (1975)  DMSM  Project  - High  Latitude  Aircraft  HF 
Propagation  Experiment,  AFCkL-TR-7^-0290,  Environmental  lies.  Papers,  No.  516 
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TABLE  12 


Digicoder  Integers  Coherence  Values 


0 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 


Coh  < 0.3095 
0.3095  £ Coh  < 0.5649 
0.5649  < Coh  < 0.7443 
0.7443  < Coh  < 0.8565 
0.8565  < Coh  < 0.9217 
0.9217  < Coh  < 0.9579 
0.9579  < Coh  < 0.9776 
0.9776  < Coh  < 0.9881 
0.9881  < Coh  < 0.9937 
0.9937  < Coh  < 0.9967 
0.9967  < Coh  < 0.9982 
0.9982  < Coh  < 0.9991 
0.9991  < Coh  < 0.9995 
0.9995  < Coh  < 0.9997 
0.9997  < Coh  < 0.9999 
0.9999  < Coh 


; 
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